Method and apparatus for tiling light sheet selective plane illumination microscopy with real-time optimized light sheet

ABSTRACT

Disclosed are a method and apparatus for tiling light sheet selective plane illumination microscopy (TLS-SPIM) with real-time light sheet optimization. The method and apparatus image multi-cellular specimens in 3D, with an improved 3D imaging ability of SPIM in resolving complex structures and optimizes SPIM live imaging performance by using a real-time adjustable tiling light sheet and creating a flexible compromise between spatial and temporal resolution. A 3D live imaging ability is provided of the TLS-SPIM by real-time imaging cellular and sub-cellular behaviors.

PRIORITY

This application claims priority to Provisional Patent Applications No. 62/119,552, 62/150,531, 62/198,969 and 62/220,572 filed Feb. 23, 2015, Apr. 21, 2015, Jul. 30, 2015 and Sep. 18, 2015, respectively, the content of each of which is incorporated herein by reference.

BACKGROUND OF THE INVENTION

This invention relates generally a fluorescence microscopy technique to three dimensional (3D) imaging of live biological specimens in real time.

Cells live in 3D environments. A more accurate understanding of cell behaviors and cell-cell interactions can be obtained by studying cells in their native, multi-cellular environments than when cultured on substrates. In a multi-cellular organism, the intracellular activities of a cell affect not only the cell itself but also its neighbor cells, and the cellular behavior of a cell results from both its own intracellular activities and its interactions with other cells. In other words, cell behaviors in a multi-cellular process are caused by the intracellular activities of all cells involved in the process. Thus, in order to understand cell behaviors in a multi-cellular process, it is necessary to study all involved cells at sub-cellular level to acquire the dynamic information of their intracellular activities, by which the underlying connections between cellular behaviors and intracellular activities of the involved cells can be revealed. For this reason, fluorescence imaging techniques that allow low-invasive 3D imaging of multi-cellular specimens with high spatial and temporal resolution are required.

Multi-cellular specimens are difficult to image in 3D because of their size and complexity accompanied with optical aberrations and light scattering. Both high 3D spatial resolution and good optical sectioning capability are required in a large field of view (FOV) in order to visualize the specimen with sub-cellular, or even cellular structural details. Meanwhile, the photobleaching and photodamage must be low enough to allow live imaging at needed speed for a certain period of time. Therefore, selective plane illumination microscopy (SPIM) is getting increasing attention for its advanced 3D live imaging ability [1-4]. Generally, the latest SPIM techniques allow low-invasive, high-speed 3D imaging of either single cell specimens with sub-cellular level, submicron 3D spatial resolution [5-8], or multi-cellular specimens with cellular level spatial resolution of a few microns [9-12], although the actual performance varies depending on the sample and fluorophore. However, imaging multi-cellular specimens in 3D with sub-cellular spatial resolution remains challenging despite the progress that has been made in SPIM development.

SPIM obtains 3D imaging ability by confining the illumination light near the detection focal plane with a light sheet. With a given detection numerical aperture (NA), its 3D imaging ability, including the spatial resolution, optical sectioning capability, field of view (FOV), photobleaching and photodamage, is mainly determined by the intensity profile of the light sheet [13, 14]. A uniformly thin and large light sheet is therefore required in SPIM to maximize its imaging ability, and the generation of such a light sheet has been a major focus of SPIM development [5-7, 15, 16]. Although different methods have been developed, none of them are ideal. Essentially, every light sheet balances the properties of light sheet thickness, the illumination light confinement, and the light sheet size differently, which results in different 3D spatial resolution, optical sectioning capability and FOV, respectively [13, 14]. Nevertheless, it becomes extremely difficult to balance these properties as the FOV increases to image multi-cellular specimens of dozens of microns or larger, since an ideal light sheet that has thin thickness, good light confinement and a large size at the same time does not exist due to the diffraction of light. Either the spatial resolution, optical sectioning capability, or both must be sacrificed to reach a larger FOV because the light sheet either becomes thicker, or the excitation light is less confined as its size increases. Therefore, the tradeoff between the spatial resolution, optical sectioning capability and FOV sets a fundamental limit on conventional SPIM, and a key problem of imaging multi-cellular specimens with sub-cellular spatial resolution using SPIM turns out to be how to increase the FOV without losing the spatial resolution and optical sectioning capability. Different approaches other than finding a perfect light sheet must be made.

Multiview SPIM is a different approach that works in two ways to improve the 3D imaging ability of SPIM on multi-cellular specimens. First, by sending the excitation light sheet and collecting the fluorescence signal from different directions, the final image is less affected by the optical aberration and light scattering introduced by the sample [9-12]. Next, 3D images taken from different directions can be fused together with similar methods used in tomography, by which the 3D spatial resolution can exceed the resolution limit set up by the light sheet in theory [1, 8, 12]. Nevertheless, the improvement in resolution by this approach relies on the number of different view directions, and the spatial resolution and the signal to noise ratio (SNR) of the 3D images taken in these different views. In multiview SPIM, the number of view directions is limited to two (lateral and axial) without rotating the sample, while the spatial resolution and SNR of the 3D image taken in each view are still strictly constrained by the light sheet used to take the image. Therefore, multiview SPIM cannot bypass the fundamental tradeoff of SPIM set up by the light sheet. Furthermore, high-speed sub-cellular dynamics, the optical aberration and light scattering introduced by both the sample and the agarose gel usually used to mount the sample make it more difficult to fuse different views accurately with high spatial resolution. Consequently, multiview SPIM is suited for imaging relatively large and bright multi-cellular specimen with cellular level spatial resolution. Imaging multi-cellular specimens with sub-cellular spatial resolution remains a problem with multiview SPIM.

The nonexistence of a perfect light sheet not only limits the 3D imaging ability of SPIM, but also creates a severe yet underappreciated practical problem. Optimization of the light sheet, including its type and dimensions, is required case by case in SPIM depending on the sample to be imaged and the biological question to be studied, as different light sheets balance the 3D imaging ability of SPIM differently [14]. The light sheet is usually determined based on the prior knowledge of the specimen and the biological process to be imaged, and it remains the same during the entire imaging process because a realignment of the microscope, that often takes hours even for SPIM experts, is usually required to change the light sheet. It is not only extremely inconvenient, but also prevents SPIM from adapting itself to reach the optimal imaging performance, because different organelles have different structures and dynamics, and both of these can change in live imaging. Instead, the optimization of SPIM imaging performance would become a closed-loop process if the light sheet could be adjusted in real-time using the immediate imaging result as feedback. More importantly, the understanding of a biological process relies on the spatial and temporal information that can be obtained from the process, whereas an ideal method that is capable of acquiring both with a high resolution is often unavailable, especially for multi-cellular specimens. The only option under such circumstances is to compromise either the spatial or the temporal resolution strategically, and only acquire the necessary information needed to understand a biological process. Such ability and flexibility to compromise can only be obtained by adjusting the light sheet in real-time in SPIM. Thus, be able to optimize the excitation light sheet in real-time is another key to improving the 3D live imaging ability of SPIM on multi-cellular specimens.

SUMMARY OF THE INVENTION

To overcome shortcoming of conventional SPIM in imaging large multicellular specimens, provided is a method of three-dimensional (3D) multi-cellular specimen imaging, with the method including using a real-time adjusted, tiling light sheet (TLS) selective plane illumination microcopy, wherein both high spatial resolution and good optical sectioning capability are maintained within a large field of view (FOV), and wherein the FOV is larger than the tiling light sheet.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

The following detailed description of preferred embodiments is made with reference to the accompanying drawings. In describing the invention, explanation of related functions or constructions known in the art is omitted for the sake of clarity in understanding the concept of the invention and to avoid obscuring the description of the invention with unnecessary details.

Recently, we described a tiling light sheet SPIM technique (TLS-SPIM) to overcome the SPIM tradeoff between the spatial resolution, optical sectioning capability and FOV by tiling the excitation light sheet [17]. Instead of using a large light sheet to illuminate the entire FOV simultaneously, a small but thin light sheet is tiled quickly within the image plane to illuminate the whole FOV, while only the fluorescence signal generated at center of each tiled light sheet is used in the final image. As a result, both high spatial resolution and good optical sectioning capability are maintained within a FOV that is much larger than the light sheet itself. Potentially, multi-cellular specimens could be imaged by TLS-SPIM with either higher spatial resolution, better optical sectioning capability or both compared to regular SPIM. Besides TLS-SPIM, several other methods were also developed to increase the FOV of SPIM with different approaches [18-20]. Unfortunately, the imaging ability of these methods are limited by using either two-photon excitation or optical components with slow dynamic performance, while none of these methods can implement all of the latest SPIM light sheets and optimize the light sheet in real-time.

Convinced by the imaging ability of TLS-SPIM observed in our previous research [17], we continued to develop and explore the potential of this method to facilitate cell behavior study in multi-cellular specimens. We present TLS-SPIM with real-time light sheet optimization to address the above problems that limit SPIM in imaging multi-cellular specimens with sub-cellular spatial resolution. We extended the 3D live imaging ability of TLS-SPIM by enabling the implementation, tiling and real-time optimization of the latest SPIM light sheets in TLS-SPIM, including the Gaussian, Bessel and Lattice light sheet. We compared the 3D imaging ability of TLS-SPIM in resolving complex structures with regular SPIM methods and verified its 3D live imaging performance on different multi-cellular specimens. In addition, we introduce a methodology to study cell behaviors in multi-cellular specimens using TLS-SPIM, which is to image intracellular activities of the cells involved in a multi-cellular process, identify their intrinsic connections with the cellular behaviors, and determine their impact on the multi-cellular process. We demonstrate the methodology using C. elegans embryo left-right symmetry breaking behavior as an example.

Our new TLS-SPIM microscope was designed to satisfy two requirements. First, the microscope can implement the latest SPIM light sheets to use their advantages in different applications, and the light sheet can be adjusted quickly to optimize the TLS-SPIM imaging performance in real-time. Next, the light sheet can be tiled rapidly within the detection image plane to enlarge the FOV without affecting the spatial resolution and optical sectioning capability. To achieve these goals, two binary spatial light modulators (SLM) are implemented in sequence in our TLS-SPIM microscope to modulate the excitation light. The two SLMs are conjugated to the image focal plane and the rear pupil of the excitation objective, respectively. Either a Gaussian light sheet, a Bessel light sheet or a Lattice light sheet can be used in the microscope, and it takes less than a millisecond to either tile or change the light sheet by applying different binary phase maps to the SLMs.

To examine the ability of TLS-SPIM in resolving complex structures, we imaged a C. elegans embryo (OD95) expressing GFP::PLC∂PH (membrane) and H2B::mCherry (nucleus) every 30 minutes for ˜5 hours, from ˜50 cell stage to comma stage until the muscle cells were functioning. A ˜0.7 μm thick, ˜10 μm long Bessel light sheet, generated by scanning a single Bessel beam, was tiled at three positions at 10 μm intervals to image the embryo. The nearly isotropic spatial resolution (˜320 nm lateral, ˜460 nm axial) and good optical sectioning capability enabled the identification of individual cells at most developmental stages. Such imaging ability will allow the study of morphogenesis at the necessary sub-cellular resolution through later stages of embryonic development than was previously possible via other techniques. Specifically, we show that it is difficult to obtain similar results by conventional SPIM with either Gaussian, Bessel or Lattice light sheets, as the light sheet must be thicker or the illumination light must be less confined to image the same FOV in regular SPIM. The comparison also shows that the imaging performance of different light sheets varies from sample to sample due to the different sample properties. Instead of sticking to a particular light sheet or method, a more suitable light sheet, including the type, dimensions, tiling number and tiling positions, can always be implemented in TLS-SPIM in real-time with the immediate feedback from the imaging result.

To demonstrate the 3D imaging ability of TLS-SPIM on larger specimens, we imaged the tailbud of a ˜15 hpf nuclei-labeled zebrafish embryo in the next. A ˜1.2 μm thick, ˜30 μm long Bessel light sheet, which provides ˜320 nm lateral resolution, ˜750 nm axial resolution, was tiled at nine positions at ˜25 μm intervals to image a 0.2 mm×0.2 mm volume. Most cell nuclei can be identified and examined with sub-cellular structural details. Again, the structure was better resolved by TLS-SPIM as compared to regular SPIM attributed to the usage of a thinner light sheet with better light confinement. Nevertheless, as can be observed from the result, the 3D imaging ability of TLS-SPIM on multi-cellular specimens is still limited by optical aberration and light scattering that affect all light microscopy techniques.

Due to the light sheet tiling in TLS-SPIM, each image plane is illuminated multiple times equivalent to the number of tiling positions, which raises concerns of slower imaging speed and higher photobleaching and photodamage. In regular SPIM, either the spatial resolution, optical sectioning capability or both must be traded for a larger FOV, while in TLS-SPIM, the temporal resolution can be traded as a substitute for FOV by tiling the light sheet. The decision of whether to sacrifice the temporal resolution, and how much to sacrifice will depend on the biological question to be answered. In addition, the problems caused by light sheet tiling are fully predictable and controllable by limiting the number of tiling positions.

To evaluate the live imaging performance of TLS-SPIM, we imaged endogenously labeled myosin II particle activities in a live C. elegans embryo (LP162) every 8.5 seconds for 150 time points. A Bessel light sheet, generated by scanning an incoherent seven Bessel beam array, was tiled at three positions to image the embryo. Both the spatial and temporal resolution of our results are sufficient to track and analyze the myosin II particle movements in 3D, by which the behavior of the actomyosin network and its impact on embryonic cells can be studied. In another example, we imaged a C. elegans embryo (OD95) with the same light sheet every 30 seconds in two channels for 167 time points. Both high spatial resolution and good SNR were maintained through the process, and the embryo hatched successfully in the end. The result shows comparable spatial resolution, SNR, imaging speed and much longer imaging time compared to the results obtained from similar specimens imaged by Bessel SPIM and Lattice light sheet microscopy in previous publications [6, 7]. The result allows us to extract the dynamic information of every individual cell and analyze their relationship over developmental time, facilitating both visualization and understanding of cell behaviors as embryonic development progresses. Clearly, TLS-SPIM performs well for 3D live imaging in terms of the spatial resolution, SNR, imaging speed, photobleaching and photodamage when the tiling number is kept low.

A new methodology can be established to study cellular behaviors in multi-cellular specimens with the dynamic information of intracellular activities acquired by TLS-SPIM. One area of active interest has been the identification of the left-right symmetry breaking mechanism of C. elegans embryos at the 4 to 6 cell stage [21-24]. TLS-SPIM can provide high spatial and temporal resolution information of different intracellular activities involved in the process, such as the actomyosin network behavior and the deformation of each individual cell, implicating the mechanical mechanism that drives symmetry breaking. We identified three intracellular activities that could have major contributions to the symmetry breaking behavior mechanically. First, the torque generated by the unaligned contractile rings of the ABa and ABp cells during mitosis could shear the ABar and ABpr to ABal and ABpl cells. Second, the simultaneous relative rotation of their daughter cells, ABal to ABar and ABpl to ABpr, could be caused by the torque generated by the counter flow of asymmetrically distributed myosin II particles on the contractile ring of ABa cell, which is consistent with previous observation [23]. Third, the retreating of the EMS cell, entering mitosis afterwards, could assist the symmetry breaking behavior and the positioning of all daughter cells by pulling the neighbor cells. Our observations not only provide evidence to verify some of the previous hypothesis [23, 24], but also suggest new mechanism to explain the behavior. Based on these observations, more specific experiments can be designed, quantitative image analysis tools and mechanical modeling methods can be introduced to analyze and understand the biochemical and biophysical mechanism of the symmetry breaking behavior. As demonstrated, with the advanced 3D live imaging ability of TLS-SPIM, instead of treating different cell behaviors isolatedly, a methodology of imaging and analyzing intracellular activities of different cells and finding their intrinsic connections can be applied to study cell behaviors in multi-cellular specimens.

The methodology can certainly be applied to later stage C. elegans embryos and other multi-cellular specimens with the 3D imaging ability of TLS-SPIM. For instance, by examining all cells of the same embryo at two time points of 1.5 minute apart at the 8 to 12 cell stage, we noticed that both the shape and position of the C cell, E cell and MS cell changed following the mitosis of the ABal, ABar, ABpl and ABpr cells, that divided almost synchronously. The positions of ABal, ABar, ABpl and ABpr daughter cells also seem affected by the synchronized mitotic processes of these cells. The consistent observations of aggressive cell position rearrangement during synchronized mitotic processes through the embryo development also suggest that the contractile rings of synchronized mitotic cells play an important role in cell positioning and morphogenesis of early stage C. elegans embryos. More imaging and quantitative analysis on different intracellular activities that are related to force generation are required to understand these behaviors completely.

TLS-SPIM requires more tiling positions to further increase the FOV and maintain the same spatial resolution and optical sectioning capability, making it less ideal for high resolution live imaging of large multi-cellular specimens. As no imaging technique can achieve submicron 3D spatial resolution and high imaging speed simultaneously on specimens of hundreds of microns or larger, we developed a compromised solution to handle this problem, which is to acquire the spatial and temporal information of a biological process separately, namely, alternating the 3D imaging between a high spatial resolution, low imaging speed mode and a low spatial resolution but high imaging speed mode. In this way, the results acquired in different modes could still provide sufficient spatial and temporal information to understand the imaged process, as long as the switching between different modes is fast enough. Our TLS-SPIM microscope enables a quick switch by adjusting the excitation light sheet in less than a millisecond, realized by changing the phase maps applied to the binary SLMs used in the microscope. This aspect was demonstrated by imaging cell migration of mosaically labeled mesodermal cells expressing a GFP-tagged membrane marker in the tailbud of a ˜15 hpf zebrafish embryo. The region of interest (80 μm×80 μm×50 μm) was imaged with ˜320 nm lateral resolution, ˜750 nm axial resolution every ˜9 seconds, and with ˜320 nm lateral resolution, ˜450 nm axial resolution every minute alternately. By this method, both the 3D structure and the dynamics of the migrating cells can be clearly visualized. A notable difference between the imaged cells and cells migrating on 2D substrates is the absence of membrane ruffling behavior, despite the similar lamellar protrusions [7], which could be caused by more strict constraints on cell membrane by neighbor cells in multi-cellular specimens. A reasonable understanding of a biological process can still be obtained by compromising the imaging ability flexibly when an ideal method is not available. Such ability is especially valuable in imaging large multi-cellular specimens as they are more difficult to image and the cellular behaviors are more complicated in such specimens.

TLS-SPIM with real-time light sheet optimization not only offers advanced 3D live imaging ability, allows new methodology to be applied to understand cell behaviors in multi-cellular specimens, but also enables the delivery of such imaging ability and analysis method to general SPIM users with its special instrument design, as there is no hardware change required to operate the microscope after the initial alignment, including the light sheet optimization, tiling, and imaging mode switching. A deep understanding on TLS-SPIM technique itself is not necessary for general users to use this technique, the most suitable imaging condition, including the light sheet and its dimensions, tiling numbers, positions and imaging modes can be determined solely based on the sample features, desired spatial resolution, SNR, imaging speed, FOV and the acquired images. The phase maps to be applied to the SLMs can be calculated automatically and applied accordingly. On the other hand, TLS-SPIM can take the advantages of most previously developed SPIM techniques rather than conflict with any of them, such as the different light sheets and the multiview SPIM configuration. Altogether, TLS-SPIM with real-time light sheet optimization will make SPIM much more feasible and favorable than it used to be to general biologists.

Two perceived problems are still limiting the 3D imaging ability and application of SPIM. First, SPIM is still limited by optical aberration and light scatting as other optical imaging techniques, and the problem is doubled in SPIM due to the separated illumination and detection. An effective combination of SPIM and adaptive optics without heavily affecting its live imaging performance is desired. After all, how well SPIM works relies on how well the excitation light is confined. Next, SPIM is unique due to the enormous spatial and temporal information that can be collected with it. However, appropriate quantitative image analysis tools are not yet widely available to general users to make most of the collected information useful. On the other hand, the development of appropriate image analysis tools itself faces many challenges at the same time, because of the huge variety of the biological questions, the imaging ability difference of the different SPIM techniques, and the lack of clear indications of what kind of quantitative information is desired and how the information can be used. The gap between biologists, physicists, and computer scientists must be filled to have this problem solved. Ultimately, moving cell biology research from 2D to 3D is inspiring but still challenging. The solutions to the above problems will perhaps determine how much impact SPIM can make in accelerating the transition.

Methods Design Considerations

The Key Tradeoff in SPIM:

Axial resolution, optical sectioning capability and FOV are the three key factors affecting the 3D imaging ability of SPIM. These factors are only determined by the intensity profile of the excitation light sheet in regular SPIM with a given detection NA. The axial resolution, defining the smallest axial separation SPIM can resolve in theory, relies on the thinnest axial component of the excitation light sheet in real space, which corresponds to the highest axial frequency of the light sheet OTF in frequency space. The optical sectioning capability measures the ability of SPIM to suppress out-of-focus fluorescence background. It determines whether the theoretical spatial resolution, especially the axial resolution, can be achieved with high enough SNR in practice. The optical sectioning capability of SPIM depends on the percentage of the illumination light confined within the detection depth of focus in real space, and it is also reflected by how quickly the light sheet OTF amplitude rolls off to zero along the axial direction in frequency space [13, 14].

Aside from the excellent axial resolution and optical sectioning capability of SPIM, FOV is a major limiting factor of SPIM due to the optical configuration of SPIM. In general, it is limited by the length of the excitation light sheet along the light propagation direction, where the light sheet intensity profile remains uniform. For this reason, the light sheet length must be increased to enlarge the SPIM FOV for large specimens. However, no matter what kind of light sheet is used, the light sheet either becomes thicker or the illumination light is less confined as its length increases due to the diffraction of light. As a result, the size of SPIM FOV always goes opposite to the SPIM axial resolution and optical sectioning capability. Therefore, the relationship between SPIM axial resolution, optical sectioning capability and FOV represents a key tradeoff in SPIM, and this tradeoff cannot be overcome by regular SPIM due to the diffraction of light. TLS-SPIM provides a solution to this problem.

Real-Time Optimization of the Excitation Light Sheet:

Because of the tradeoff in SPIM, no light sheet is suited for all specimens. Optimization of the light sheet is required in SPIM, and the benefit of SPIM is not guaranteed unless an appropriate light sheet is selected based on the specimen to be imaged and the question to be answered. Different light sheets have been developed for SPIM imaging. Typical light sheets include Gaussian light sheets, Bessel light sheets and Lattice light sheets, etc [1-7, 15]. For light sheets of comparable sizes, Gaussian light sheets have the most restricted light confinement, but they are usually too thick to obtain submicron axial resolution. Bessel light sheets contain a much thinner central peak that enables high axial resolution to be obtained in theory, but the optical sectioning capability is sacrificed due to the poor illumination light confinement. Lattice light sheets reach a better balance between light sheet thickness and light confinement compared to the previous two by using a coherent Bessel beam array, but it is more sensitive to optical aberration and light scattering. [13, 14]

On the other hand, SPIM lowers the photobleaching and photodamage compared to conventional microscopy techniques. Obviously, light sheets confining more excitation light within the detection depth of focus produce less photobleaching and photodamage. Nevertheless, the photodamage of SPIM not only depends on the light sheet intensity profile, but also on how it is generated. In SPIM, the excitation light sheet can be either a real one, existing simultaneously across the entire illumination field, or a virtual one, created by scanning an excitation beam or dithering a beam array. As a consequence, for light sheets of the same intensity profile, a real light sheet has a lower peak intensity than a virtual one created by dithering a beam array or scanning a single excitation beam. Although it hasn't been studied systematically due to the complexity of a fair comparison and different natures of different specimens, lower excitation peak intensity is usually preferred to decrease the photodamage. Using a light sheet with lower peak intensity should always be considered when photodamage becomes a major problem [6, 7]. A general guidance of how to choose SPIM light sheet can find in previous publications [13, 14]. For the above reasons, to adjust the light sheet in real-time based on the immediate feedback from the imaging result is important to improve the 3D live imaging ability of SPIM. We conceived our research and constructed our microscope accordingly.

Microscope Configuration

CW lasers with excitation wavelengths of 488 nm and 561 nm are used for linear excitation (Coherent, Sapphire LP, 200 mW). Laser beams from both lasers are expanded to a common 1/e² beam diameter of 1.5 mm and combined into a single co-linear beam using a LaserMux dichroic beam combiner. The combined laser beam is sent to an acousto-optical tunable filter (AOTF, AA Opto-Electronic, nAOTFnC-400.650-TN) used to select one or more wavelengths and modulate the laser beam intensity. The laser beam is directed to one of the two separated optical paths after the AOTF and a half wave plate HWP1.

The microscope is operated in two modes by sending the laser beam to the two separated optical paths. In the first mode, a binary spatial light modulator (SLM), SLM2, conjugated to the rear pupil of the excitation objective is used to generate a single excitation beam or an incoherent beam array. A virtual excitation light sheet is generated by scanning the excitation beam or dithering the excitation beam array using a mirror galvanometer. SLM2 is also used to tile the excitation light sheet by changing the excitation beam/beam array position along the beam direction. In the second mode, Lattice light sheets are used for SPIM imaging. Two binary SLMs, SLM1 and SLM2 are implemented in sequence, in which SLM1 is conjugated to the image focal plane of the excitation objective and SLM2 is still conjugated to the rear pupil of the excitation objective. A Lattice light sheet is generated by applying the corresponding phase maps to both SLM1 and SLM2, and the light sheet was tiled by SLM2 in the same way as that in the first mode.

In the second mode, the laser beam is expanded to 15 mm in vertical direction by a pair of cylindrical lenses (focal length CL1=25 mm, CL2=250 mm) and sent to SLM1. The SLM assembly consisted of a polarizing beam splitter cube, a half-wave plate and a 1280×1024 pixel binary SLM (Forth Dimension, SXGA-3DM). The SLM1, conjugated to the excitation objective image focal plane, is imaged to the second identical binary SLM assembly, SLM2, through lens L5=500 mm. After that, the laser beam is sent to the beam scanning assembly consisted of a pair of galvos (Cambridge Technology, 6215HP-1HB) for laser beam and light sheet scanning. Both galvos are conjugated to SLM2, through relay lenses L6=150 mm, L7=75 mm, and conjugated with each other through relay lenses L8=100 mm, L9=100 mm, and finally conjugated to the rear pupil of the excitation objective (Nikon, CFI Apo 40×W NIR 0.8 NA) through relay lenses L10=75 mm and L11=150 mm. An optical slit was placed at the focal plane of lens L6, to block the undesired diffraction orders generated by SLM2. The detection objective (Nikon, CFI Apo 40×W NIR 0.8NA), mounted on an objective scan piezo (PI, P-724), is placed orthogonal to excitation objective. The emitted fluorescence is collected through the detection objective and imaged onto a sCMOS camera (Hamamatsu, Orca Flash 4.0) by tube lens L12=225 mm. In the first operation mode, the laser beam is directed to the other optical path by flip mirrors FM1 and FM2. The laser beam is expanded to 6 mm diameter by relay lenses L3=25 mm and L4=100 mm and sent to SLM2 after that. The laser beam followed the same optical path as that in the other mode after SLM2. Manual operation is required to switch the microscope between the two modes. The flip mirrors FM1 and FM2 can be replaced by galvos to allow fast switching between the two operation modes. A wildfield optical path used to find samples before SPIM imaging is not shown in the schematic diagram. The control system and hardware are the same as what reported in the previous publications [5, 6, 13]. The sCMOS camera triggers are used to synchronize the operation of the SLMs, galvos, objective piezo and camera itself.

Microscope Working Principle

The microscope was designed to satisfy two requirements. First, the excitation light sheet can be tiled quickly within the detection image plane to enlarge the FOV without affecting the axial resolution and the optical sectioning capability. Next, both the type and the dimension of the light sheet can be changed quickly to optimize the SPIM imaging performance in real-time.

Tiling of the excitation light sheet can be realized by adding a variable spherical phase map to the excitation light wavefront at a plane conjugated to the excitation objective rear pupil. It can be achieved by using either a SLM, a deformable mirror, or a focus variable lens. We use a binary SLM because it is flexible, easy to control and the light sheet can be positioned fast and accurately by applying different phase maps. The applied phase map can be refreshed in less than a millisecond (˜3.2K Hz rate), which enables both fast tiling and real-time optimization of the light sheet. The method of tiling the light sheet with a binary SLM by applying a spherical phase map was described in details in our previous publication [17].

The excitation light sheet used in SPIM is usually generated by scanning a single excitation beam or dithering a beam array. The intensity profile of light sheet is adjusted by changing the beam or the beam array used to create the light sheet. It can be realized by modifying the intensity and phase of the excitation laser beam. For example, circular laser beams of different diameters, corresponding to different excitation numerical apertures (NA_(exc)) can be used to create Gaussian beams of different dimensions. Annular laser beams of different inner and outer diameters, corresponding to different inner and outer excitation NA (NA_(ID), NA_(OD)), can be used to create different Bessel beams. A metal coated quartz mask containing transmissive patterns of different geometries is usually used to shape the laser beam intensity profile and create these beams [5-7, 13, 14]. However, the number of options is limited and realignment is required to switch between different beams in such a way.

In our microscope, SLM2 is also used to modify the intensity and phase profile of the excitation laser beam to adjust the light sheet intensity profile in real-time besides tile the light sheet. The SLM2 is separated to two areas, retain area and abandon area. A high frequency binary phase grating is applied to the abandon area, so that the incident light falling on the abandon area is diffracted to different directions from the light falling on the retain area. The abandoned light is lost or blocked by the optical slit afterwards. The retained light is allowed to propagate along the designed optical path. For example, the annular amplitude mask used to generate a Bessel beam can be replaced by applying the phase map to SLM2. It is calculated by multiplying a reverse binary amplitude mask, and a high frequency binary phase grating (2 pixels per period). In this way, the laser beam intensity profile can be modified by adjusting the shape and size of the retain area. Further modification can be done by controlling the phase within the retain area of SLM2. Most excitation beams/light sheets can be generated by this method, and we developed specific solutions to generate either a single Gaussian/Bessel beam, an incoherent Gaussian/Bessel beam array or a coherent Bessel beam array in our microscope, so that both Gaussian light sheets, Bessel light sheets and Lattice light sheets can be used in our microscope. The details are discussed below.

Single Excitation Beam

The microscope is operated in the first mode, and only SLM2 is used. The excitation beam and its dimension are controlled by shape and size of the retain area on SLM2. The tiling of the light sheet is realized by applying a binary spherical phase map to the retain area. Extra diffraction orders are blocked by the optical slit. The details to generate the binary spherical phase map is discussed in our previous publication [17]. For example, a Bessel beam (NA_(OD)=0.35, NA_(ID)=0.14) was generated by applying the phase map. The same beam was moved to the right side of the previous position by applying the phase map. A thicker and longer Bessel beam (NA_(OD)=0.2, NA_(ID)=0.05) was generated on the left side of the initial beam position by applying the phase map.

Incoherent Beam Array

The microscope also works in the first mode to generate an incoherent excitation beam array. A light sheet with lower peak intensity can be generated by dithering the beam array. An incoherent beam array with beams of roughly equal intensity can be generated by using a Dammann grating [25, 26], which is a periodic binary phase pattern. In the meantime, the beam array can be tiled by superimposing a continuous spherical phase map to the Dammam grating. A binary phase map to generate the beam array at a desired position is obtained by resetting the pixel values of the superimposed phase map to 0 and π. Values between 0 and π were set to 0, and values between π and 2π were set to π. The final binary phase map applied to SLM2 is obtained by replacing the phase map of the abandon area with the high frequency grating pattern.

However, this method has several limitations due to the limited SLM pixels. To create a beam array of N beams using a Dammann grating, 2×int[(N−1)/4]+2 transition points are required per grating period [26], while the minimum distance between any two adjacent transition points must be larger than one pixel to allow the grating pattern to be applied to the SLM. In consequence, the pixel number per grating period has to be increased as the number of beams increases, which decreases the corresponding beam array period at the same time. On the other hand, the beam array period has to be large enough to make the beam array wide enough to cover the desired FOV and avoid the interference between different beams, because the Dammann grating doesn't control the phase relationship between different beams. For above reasons, it is difficult to generate a beam array of more than 21 beams by this method, and the selection of beam array period corresponding to each beam number is also limited.

Coherent Bessel Beam Array

TLS-SPIM works in the second mode to use a tiling Lattice light sheet. Both SLM1 and SLM2 are used to generate and tile a coherent Bessel beam array in this mode. As reported in previous research from Betzig's lab [7], a binary phase map generated from the desired beam array amplitude and phase profile is applied to a binary SLM that conjugated to the image focal plane of the excitation objective. The SLM is imaged to a transmissive annular pattern on a photo mask, which is conjugated to the excitation objective rear pupil to remove the undesired diffraction orders. A similar method is used in our microscope to generate a coherent Bessel beam array except that SLM2 is used to replace the photo mask and tile the Lattice light sheet. A binary phase map is applied to SLM1 to generate a coherent excitation beam array. The high frequency grating pattern is also used to remove the undesired illumination light. Two phase maps applied to SLM2 (NA_(OD)=0.35, NA_(ID)=0.14) to shape the excitation laser beam and tile the coherent Bessel beam array.

A beam array that is more close to an incoherent beam array can also be generated in this mode with the same method by increasing the beam array period. The beam array period and the number of beams can be controlled with more flexibility in this mode compared to the previous one, despite the more complicated optical system and more laser power loss due to the use of two SLMs.

Microscope Alignment and Calibration

The microscope alignment is similar to what reported in our previous publication [13]. The binary phase maps used to generate different excitation beams/beam arrays at different positions in the two working modes were calculated in MATLAB in advance and calibrated in fluorescent dye solution. All phase maps were loaded into the corresponding SLM prior to experiments. Selected phase maps are loaded in sequence and synchronized by the sCMOS camera during 3D imaging. Potentially, phase maps corresponding to different light sheets can be calculated in-situ with the given parameters of the optical system components after the calibration.

Image Analysis

3D image stacks collected at different light sheet tiling positions were stitched together according to the length and tiling positions of the light sheet. The combined image stack was deconvolved using the Richardson-Lucy 3D deconvolution algorithm in MATLAB [13]. The point spread function (PSF) corresponding to the light sheet was collected by imaging an isolated 100 nm fluorescent particle. The deconvolved result was filtered with low pass filters to remove noises beyond the microscope OTF. Cell segmentation was performed using Matlab and Amira. Analyzed image data were rendered in 3D with Amira for display. More accurate image registration method to combine different image stacks is under development.

Sample Preparation

C. elegans embryos were dissected from adult worms. Dissected C. elegans embryos were placed directly on top of poly-D-lysine-coated coverslips for imaging. The following strains were used: OD95 1tIs37 [pAA64; pie-1::mCherry::HIS-58+unc-119(+)] IV. 1tIs38 [pAA1; pie-1:GFP::PLC∂^(PH)+unc-119(+)], LP162 nmy-2(cp13[nmy-2::gfp+LoxP]) I. RW10029 zuIs178 [his-72(1 kb 5′ UTR)::his-72::SRPVAT::GFP::his-72 (1 KB 3′ UTR)+5.7 kb XbaI−HindIII unc-119(+)]. stIs10024 [pie-1::H2B::GFP::pie-1 3′ UTR+unc-119(+)].

Nuclei of zebrafish embryos were labeled by ubiquitous expression of nuclear localized Kikume green/red photoconvertible protein (NLS-KikGR). Membrane labeling was achieved by injecting 25 picograms of a reporter plasmid at the 1-cell stage. The plasmid consisted of a 1 kb fragment of the ntla promoter driving expression of membrane localized Cherry, in order to enrich for expression in mesodermal cells. Injected plasmid DNA is inherited mosaically in zebrafish embryos, causing cells to be scatter-labeled. Embryos were imaged in their chorions at the 6-12 somite stages, focusing on the posteriorly localized embryonic structure of the tailbud. To mount the zebrafish embryo, a drop of 1% low melting temperature agarose gel was dropped on top of a cleaned coverslip, and the zebrafish embryo was place on top of the melted Agarose drop afterwards. The embryo is ready for imaging after the agarose gel is solidified. All embryos were imaged in water at 20° C.

The invention is not limited to the disclosed preferred embodiments, and should be construed to cover all such alternatives, modifications and equivalents as defined in the appended claims.

REFERENCES

-   1. Huisken, J., Swoger, J., Del Bene, F., Wittbrodt, J. &     Stelzer, E. H. K. Optical sectioning deep inside live embryos by     selective plane illumination microscopy. Science 305, 1007-1009     (2004). -   2. Keller, P. J., Schmidt, A. D., Wittbrodt, J. & Stelzer, E. H. K.     Reconstruction of zebrafish early embryonic development by scanned     light sheet microscopy. Science 322, 1065-1069 (2008). -   3. Wu, Y., et al. Inverted selective plane illumination microscopy     (iSPIM) enables coupled cell identity lineaging and     neurodevelopmental imaging in Caenorhabditis elegans. Proc. Natl.     Acad. Sci. U.S.A., 108, 17708-17713 (2011). -   4. Keller, P. J. Imaging morphogenesis: technological advances and     biological insights. Science, 7, 1234168 (2013). -   5. Planchon, T. A., et al. Rapid three-dimensional isotropic imaging     of living cells using Bessel beam plane illumination. Nat. Methods     8, 417-423 (2011). -   6. Gao, L., et al. Noninvasive imaging beyond the diffraction limit     of 3D dynamics in thickly fluorescent specimens. Cell 151, 1370-1385     (2012). -   7. Chen, B. C., et al. Lattice light-sheet microscopy: imaging     molecules to embryos at high spatiotemporal resolution. Science,     346, 1257998 (2014). -   8. Wu, Y., et al. Spatially isotropic four-dimensional imaging with     dual-view plane illumination microscopy. Nat. Biotechnol. 31,     1032-1038 (2013). -   9. Tomer, R., Khairy, K., Amat, F. & Keller, P. J. Quantitative     high-speed imaging of entire developing embryos with simultaneous     multiview light-sheet microscopy. Nat. Methods 9, 755-763 (2012). -   10. Krzic, U., Gunther, S., Saunders, T. E., Streichan, S. J. &     Hufnagel, L. Multiview light-sheet microscope for rapid in toto     imaging. Nat. Methods 9, 730-733 (2012). -   11. Schmid, B., et al. High-speed panoramic light-sheet microscopy     reveals global endodermal cell dynamics. Nat. Commun. 4:2207 (2013). -   12. Chhetri, R. K., et al. Whole-animal functional and developmental     imaging with isotropic spatial resolution. Nat. Methods 12,     1171-1178 (2015). -   13. Gao, L., Chen, B. C., Shao, L. & Betzig, E. 3D live fluorescence     imaging of cellular dynamics using Bessel beam plane illumination     microscopy. Nat. Protoc. 9, 1083-1101 (2014). -   14. Gao, L. Optimization of the excitation light sheet in selective     plane illumination microscopy. Biomed. Opt. Express 6, 881-890     (2015). -   15. Fahrbach, F. O., Simon, P. & Rohrbach A. Microscopy with     self-reconstructing beams. Nat. Photon. 4, 780-785 (2010). -   16. Vettenburg, T., et al. Light-sheet microscopy using an Airy     beam. Nat. Methods 11, 541-544 (2014). -   17. Gao, L. Extend the field of view of selective plan illumination     microscopy by tiling the excitation light sheet. Opt. Express 23,     6102-6111 (2015). -   18. Zong, W., et al. Large-field high-resolution two-photon digital     scanned light-sheet microscopy. Cell Res. 25, 254-257 (2015). -   19. Dean, K. M., Roudot, P., Welf, E. S., Danuser, G. & Fiolka, R.     Deconvolution-free Sub-cellular Imaging with Axially Swept Light     Sheet Microscopy. Biophys. J., 108, 2807-2815 (2015). -   20. Chmielewski, A. K., et al. Fast imaging of live organisms with     sculpted light sheets. Sci. Rep. 5, 9385 (2015). -   21. Pohl, C. & Bao, Z. Chiral forces organize left-right patterning     in C. elegans by uncoupling midline and anteroposterior axis. Dev.     Cell 19, 402-412 (2010). -   22. Singh, D. & Pohl, C. Coupling of rotational cortical flow,     asymmetric midbody positioning, and spindle rotation mediates     dorsoventral axis formation in C. elegans. Dev. cell 28, 253-267     (2014). -   23. Naganathan, S. R., Furthauer, S., Nishikawa, M., Julicher, F. &     Grill, S. W. Active torque generation by the actomyosin cell cortex     drives left-right symmetry breaking. eLife 3, 04165 (2014). -   24. Pohl, C., Cytoskeletal Symmetry Breaking and Chirality: From     Reconstituted Systems to Animal Development. Symmetry 7, 2062-2107     (2015). -   25. Dammann, H. & Klotz, E. Coherent optical generation and     inspection of two-dimensional periodic structures. Opt. Acta 24,     505-515 (1977). -   26. Zhou, C. & Liu, L. Numerical study of Dammann array     illuminators. Appl. Opt. 34, 5961-5969 (1995). 

What is claimed is:
 1. A method of live three-dimensional (3D) multi-cellular specimen imaging, the method comprising: tiling a light sheet for selective plane illumination to enlarge a field of view (FOV) of selective plane illumination microscopy (SPIM), wherein both high spatial resolution and optical sectioning capability are maintained within the FOV, and wherein the FOV is larger than a tiling light sheet (TLS); and adjusting a dimension and an intensity of the TLS in real-time to optimize imaging performance, wherein adjusting real time imaging ability of the SPIM adjusts the TLS in real-time.
 2. The method of claim 1, wherein the adjusting of the TLS in real-time includes at least one of alternating between a high spatial resolution, a low imaging speed mode, and a low spatial resolution, high imaging speed mode.
 3. A method of real-time, three-dimensional (3D) multi-cellular specimen imaging, the method comprising: adjusting a tiling light sheet (TLS) selective plane illumination; acquiring spatial information of a biological process of the specimen; acquiring temporal information of the biological process of the specimen, wherein the spatial information is acquired separate from the temporal information; and alternating between a high spatial resolution, a low imaging speed mode, and a low spatial resolution but high imaging speed mode, wherein both high spatial resolution and good optical sectioning capability are maintained within a field of view (FOV), and wherein the FOV is larger than the TLS. 